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ABSTRACT 

The modeling of the pressure- strain correlation of turbulence is examined from a basic 
theoretical standpoint with a view toward developing improved second-order closure models. 
Invariance considerations along with elementary dynamical systems theory are used in the 
analysis of the standard hierarchy of closure models. In these commonly used models, the 
pressure-strain correlation is assumed to be a linear function of the mean velocity gradients 
with coefficients that depend algebraically on the anisotropy tensor. It is proven that for 
plane homogeneous turbulent flows the equilibrium structure of this hierarchy of models 
is encapsulated by a relatively simple model which is only quadratically nonlinear in the 
anisotropy tensor. This new quadratic model - the SSG model - is shown to outperform the 
Launder, Reece, and Rodi model (as well as more recent models that have a considerably 
more complex nonlinear structure) in a variety of homogeneous turbulent flows. However, 
some deficiencies still remain for the description of rotating turbulent shear flows that are 
intrinsic to this general hierarchy of models and, hence, cannot be overcome by the mere 
introduction of more complex nonlinearities. It is thus argued that the recent trend of 
adding substantially more complex nonlinear terms containing the anisotropy tensor may be 
of questionable value in the modeling of the pressure-strain correlation. Possible alternative 
approaches are discussed briefly. 


‘This research was supported by the National Aeronautics and Space Administration under NASA Con- 
tract No. NAS 1-1 8605 while the first and second authors were in residence at the Institute for Computer 
Applications in Science and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665. 
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1. INTRODUCTION 


The pres sure- strain correlation plays a pivotal role in determining the structure of a wide 
variety of turbulent flows. Consequently, the proper modeling of this term is essential for 
the development of second-order closure models that have reliable predictive capabilities. 
Rotta (1951) developed the first simple model for the slow pressure-strain correlation (i.e., 
the part that is independent of the mean velocity gradients) which describes the return 
to isotropy behavior of turbulence within the framework of a full Reynolds stress closure. 
This model has served as a cornerstone for the representation of the slow pressure-strain 
correlation in a variety of the commonly used second-order closures such as the Launder, 
Reece, and Rodi (1975) model. Subsequent to this work, Lumley (1978) demonstrated the 
need for nonlinear terms in models for the slow pressure-strain correlation and derived a 
nonlinear representation theorem for the slow pressure-strain correlation based on isotropic 
tensor function theory. In the high-Reynolds-number and small anisotropy limit, the Lumley 
(1978) model reduces to the Rotta model. 

The simplest model for the rapid pressure- strain correlation that is used in second-order 
closure modeling is based on the assumption of isotropy of the coefficients of the mean 
velocity gradients; this gives rise to a rapid pressure-strain model with a single term that 
is proportional to the mean rate of strain tensor (see Rotta 1972 and Mellor and Herring 
1973). Starting with the work of Launder, Reece, and Rodi (1975), anisotropic models for 
the rapid pressure-strain correlation have been formulated wherein the coefficients of the 
mean velocity gradients are taken to be functions of the anisotropy tensor. In the Launder, 
Reece, and Rodi model, the fourth-rank tensor of coefficients of the mean velocity gradient 
tensor is linear in the anisotropy tensor whereas most of the newer models developed during 
the last decade are nonlinear (see Shih and Lumley 1985, Haworth and Pope 1986, Speziale 
1987, Reynolds 1987, and Fu, Launder, and Tselepidakis 1987). The nonlinear models of 
Lumley and co-workers have been primarily developed by the use of realizability constraints 
(see Lumley 1978). In contrast to this approach, the nonlinear model of Speziale (1987) was 
derived using a geostrophic flow constraint (i.e., material frame- indifference in the limit of 
two-dimensional turbulence) whereas Reynolds (1988) has attempted to develop models that 
are consistent with Rapid Distortion Theory (RDT). 

In this paper, the general hierarchy of closure models for the pressure-strain correlation 
will be considered which are linear in the mean velocity gradients, with coefficients that are 
functions of the anisotropy tensor. This hierarchy of models, which was motivated by analy- 
ses of homogeneous turbulence, encompasses all of the closure models for the pressure-strain 
correlation that have been used in connection with second-order closures. A general repre- 
sentation for this hierarchy of closure models will be derived by means of invariant tensor 
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function theory. This general representation for the pressure-strain correlation will then be 
applied to plane homogeneous turbulent flows - the class of flows that have long played a piv- 
otal role in the screening and calibration of such models. However, there will be one notable 
difference with previous work on this subject: the simplest generic form of this hierarchy of 
models will be sought that has the same equilibrium structure in the phase space of plane 
homogeneous turbulent flows as the general model. This generic form - which will be termed 
the SSG model - is only quadratically nonlinear in the anisotropy tensor. It has the advan- 
tage of being topologically equivalent to the general model in plane homogeneous turbulence 
with the simplicity of structure that allows for the determination of all empirical constants 
based on calibrations with pertinent RDT results and two well documented physical exper- 
iments (i.e., homogeneous turbulent shear flow and the decay of isotropic turbulence). This 
new SSG model will be shown to perform better than the commonly used Launder, Reece, 
and Rodi model for a variety of homogeneous turbulent flows which include plane strain, 
rotating plane shear, and the axisymmetric expansion/contraction. However, there are still 
some remaining deficiencies in the new model, particularly for rotating shear flow. Based 
on an analysis of the bifurcation diagram for rotating shear flow, it will be shown that these 
deficiencies are intrinsic to this general hierarchy of pressure-strain models and cannot be 
eliminated by the addition of more complex nonlinear terms. The implications that these 
results have on turbulence modeling will be discussed in detail along with suggested future 
directions of research. 


2. THE GENERAL PRESSURE-STRAIN MODEL 


We will consider the turbulent flow of a viscous, incompressible fluid governed by the 
Navier-Stokes and continuity equations 


dvi dvi 

dvi 

dx{ 


dP_ 

dxi 


+ vV 2 Vi 


= 0. 


( 1 ) 


( 2 ) 


In (1) - (2), Vi is the velocity vector, P is the modified pressure, and v is the kinematic 
viscosity of the fluid. The velocity and pressure are decomposed into ensemble mean and 
fluctuating parts, respectively, as follows: 

Vi = Vi-hUi, P = P + p. (3) 

Here, the mean and fluctuating velocity are solutions of the transport equations 


dvi _ dvi 

~dt + = 


dF n dTij 

v ‘ + &: 

l/U/| l/U/j 


( 4 ) 
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( 5 ) 


dv-i _ du x dui dvi dp 2 9r,, 

~et + ^ ~ at + " v * ~ s£ 

|S-o 

OXi 

where r x j = —u x Uj is the Reynolds stress tensor. 

The Reynolds stress tensor t X] is a solution of the transport equation 

dr xj _ drij &Uj dvi dC x jk „ _ 2 

at + Vfc ax fe “ T,fc dx k Tjk !h~ k + IhT ~ Uij + £ii + " v Tij ' 


where 

C’.jfc = u x UjU k + pv~£jk + WjSik (9) 

7 a^, 5u,\ 

(10) 

_ dui du x , , 

= 2 "&rs: < u > 

are, respectively, the third-order diffusion correlation, the pressure-strain correlation, and 
the dissipation rate correlation. Equation (8) is obtained by taking the symmetric part of 
the ensemble mean of the product of the fluctuating velocity Uj with equation (6). For homo- 
geneous turbulent flows, at high Reynolds numbers where the dissipation is approximately 
isotropic, the Reynolds stress transport equation (8) simplifies to 

dvj dv x „ 2 p 

Tij ~ Tik dx k ~ Tjk ~£hT k ~ n,J + Z e£ij ( 12 ) 


where 


c = ‘"aT k aT i < 13 > 

is the scalar dissipation rate. Equation (12) becomes a closed system for the determination 
of in terms of dv x jdxj once closure models for II,-,- and e are provided. Since TUj is the 
only unknown correlation that contains directional information, it is clear that it will play 
a pivotal role in determining the structure of the Reynolds stress tensor for a given mean 
velocity field. This dominant influence of Iljj on the evolution of the Reynolds stress tensor 
in (12) has motivated researchers to rely on homogeneous turbulent flows for the testing and 
calibration of pressure-strain models. 

The fluctuating pressure p is a solution of the Poisson equation 

V 2 p_ 2 ^-— dUx ^ I dUi dUj / . . \ 

dx x dxj dxj dxi dxj dx, ' ' 
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which is obtained by subtracting the divergence of (4) from the divergence of (1). In the 
absence of boundaries, equation (14) has the general solution 


P 


= -f 

47T 7 - c 


X — X' 


8v[ , du\ du'j du' x du) 
2 ^ , + 


ifaW 

ldx\) 


(15) 


dx[ dx'j dx 1 - dx[ dx'- 

For homogeneous turbulent flows (where the mean velocity gradients are spatially uni- 
form) the pressure-strain correlation takes the form 

dv k 


Hij — Aj -(- Mij k i 


dx t 


where 


a - = — H + 

3 47 r 7-oo \9xj dxi 


dui duA dV' 


dx' k dx' t |x' 


dui _ duj 


_ , 1 f°° oui 


dui dV' 


(16) 


(17) 


(18) 


dxi ) dx k |x' — x| 

It has been shown that Aj and are functionals - over time and wavenumber space - of 

the energy spectrum tensor; see Weinstock (1981, 1982) and Reynolds (1987). In a one point 
closure, this dependence would suggest models for Aj and M^*/ that depend on the history 
of the Reynolds stress tensor and dissipation rate. The simplest such models are algebraic 
in form: 

Aj = eAi(b) (19) 

KM ijkl { b) (20) 


Mi, 


where 


tjkl 


r _ T *j 3 T kk 

b{j — ' 


'V7 


TU 

„ 1 
K = 


(21) 

( 22 ) 


are the anisotropy tensor and turbulent kinetic energy. In (19) - (20), A\j and A\{jkt can 
only depend on through b t j since they are dimensionless tensors that vanish in the limit of 
isotropic turbulence. Virtually all models for the pressure-strain correlation that have been 
used in connection with second-order closure models are of the general form of (19) - (20). 
The use of this hierarchy of models for general inhomogeneous turbulent flows is based on 
the assumption of a local homogeneous structure. 

The mean velocity gradient tensor can be decomposed into symmetric and antisymmetric 

parts as follows: 


dv x 

dxj 


— ~^ij 4 " &ij 


(23) 
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where 


(24) 


?, i = 5 {p- + P) 

2 \OXj oxi) 

i _ 

3 dxi) 


"ij = 2 


dvi 

Ox, 


(25) 


are the mean rate of strain tensor and mean vorticity tensor, respectively. Then, the model 
for the pressure- strain correlation specified by (19) - (20) can be written in the equivalent 
form 

= e4?>(b, S', S') (26) 

where 

/ K K 

% = “xj = —“a (27) 

V C 

are the dimensionless mean strain rate and vorticity tensor whereas denotes the part of 
the function fa that is linear in the mean velocity gradients and traceless. Form invariance 
under a change of coordinates requires that fa transform as 


Qf(b,S',^)Q T = f(QbQ T ,QS , Q T ,Qw'Q T ) 


(28) 


where Q is the rotation tensor (and Q T is its transpose) which characterizes a change in 
orientation of the coordinate axes. In mathematical terms, (28) requires that fa be an 
isotropic tensor function of its arguments. By using known representation theorems for 
isotropic tensor functions (see Smith 1971) to construct fa - and then by taking the linear 
and traceless part of fa - the following model for ITj is obtained: 

n,y = fa ebij + fae(b ik b k j — -b mn b mn Sij) + faK~3ij 

+faK(b ik ^j k -f bj k !>i k - -b mn ~S mn Sij) 

d (29) 

+/?5 K(bi k bidl>jt + bj k birfl>u - -bimbmn^niSij) 

+P e K(bi k uJj k + bj k uJi k ) + PrKfoihbMUji + bj k bi<juJ t i) 

where 

A = Ao (//, III) + Ai( 77 , III)-tr(b • S) + Aa( 77 , III)-tr{ b 2 • S), i= 1,2 ( 30 ) 

ft =&(//, 7 / 7 ), , = 3,4 7 

II = bijbij, III ~ bijbjkbfd 
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(31) 

(32) 


and tr(-) denotes the trace (see Appendix A for more details). Equation (29) represents the 
most general form of the hierarchy of models (19) - (20) for the pressure-strain correlation 
that is consistent with the crucial physical constraint of invariance under coordinate trans- 
formations. It will be shown later that the Launder, Reece, and Rodi model is the linear 
limit of (29) wherein fjitfhtfat and /? 6 are constants while /? 2 ,/?s, and j3 7 are zero. 

Finally, in regard to the general model, a few comments should be made concerning non- 
inertial frames of reference. In a non-inertial frame, the mean vorticity tensor must be 
replaced by the intrinsic (i.e., absolute) mean vorticity tensor defined by 

Wij = n a + e mjl n m (33) 


where O m is the rotation rate of the non-inertial frame relative to an inertial framing and 
is the permutation tensor (see Launder, Tselepidakis, and Younis 1987 and Speziale 1989a). 
Furthermore, Coriolis terms must be added to the Reynolds stress transport equation which 
then takes the form 


T *j — T ik 


&Uj 

dx^ 


^ TT , 2 X 

~ Tjk dx k n ‘ J + 3 £Sij 


+ Tjk^mki^rn) 


(34) 


in an arbitrary non-inertial reference frame. 


3. PLANE HOMOGENEOUS TURBULENCE 


We will consider the general class of plane homogeneous turbulent flows for which the 
mean velocity gradient tensor can be written in the form 


dvi 

dxj 


0 5 + ut 0 ) 

S — u 0 0 

0 0 0 / 


(35) 


Since the mean continuity equation (5) requires that dvi/dxi = —dvj j 8 x 2 in plane homoge- 
neous turbulent flows, equation (35) results by simply aligning the coordinates at a 45° angle 
relative to the principal directions of the symmetric part of dvi/dxj. Of course, in order to 
maintain the homogeneity and two-dimensionality of the mean flow, S and u> are constants 
while flj is given by 

a = (0,0,0) (36) 

(hence, the rotation is in the plane where the mean velocity gradients are applied). Equations 
(35) - (36) encompass, as special cases, plane shear, plane strain, rotating plane shear, and 
rotating plane strain turbulence. 
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The Reynolds stress transport equation (34) for plane homogeneous turbulence can be 
written in terms of the anisotropy tensor (given that t* = is the dimensionless 

time) as follows: 

(l - + 


where 


dbjj 

dt* 


e 

~sk 


"4" 2 ^ &3kj ) bj fc ( ^ik "f" ^ tk ~b ^ <j> ^3ki ) 


'V 


(0 

i 

0 \ 



0 

Uf 

S 


1 

0 

0 

v*j = 


V 

S 

0 

0 

Vo 

0 

0 / 



0 

0 

0 / 


n*; = Uij/2KS 


(37) 


(38) 

(39) 


and V = Tijdvijdxj is the turbulence production. Equation (37) must be supplemented 
with a transport model for the turbulent dissipation rate in order to obtain a closed system, 
We will consider the most standard form of the modeled dissipation rate transport equation 
given by 

„ £ OVi £ ,, M 

4 = C -'K T “d^ ~ C *K (40) 


where C7 el and C e2 can be functions of the invariants II and III of the anisotropy tensor 
(in the most commonly used form of this model, C e i and C e 2 are constants that assume the 
values of 1.44 and 1.92, respectively) see Hanjalic and Launder 1972). It will be argued later 
that some of the crucial conclusions to be drawn concerning the limitations of this hierarchy 
of closure models for the pressure-strain correlation are independent of the specific form of 
(40). Furthermore, it should be noted that virtually all of the alterations to (40) that have 
been proposed during the last decade are highly ill-behaved (see Speziale 1989b). 

The modeled dissipation rate equation (40) can be written in the dimensionless form 


Equation (41) is obtained by combining (40) with the transport equation for the turbulent 
kinetic energy 

K = V-£ (42) 


which is exact for homogeneous turbulence. When the general model for the pressure-strain 
correlation (given by (29), with aJ»j replaced by TT^j) is substituted into (37), a closed system 
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for the determination of 
solutions of the form 


b t j and ej SK is obtained. This system of equations has equilibrium 


(6,i)oc = 9ij (j, 

(43) 

( e \ (V Q\ 

(44) 

^skJoo ~ 9 \s'~s) 


where (•)<*, denotes the solution in the limit as time t — ► oo; these solutions were shown by 
Speziale and Mac Giolla Mhuiris (1989a) to attract all initial conditions. The equilibrium 
solutions (43) - (44) are obtained by solving the nonlinear algebraic equations that result 
when the time derivatives on the left-hand-sides of (37) and (41) are set to zero. It is a 
simple matter to show that there is a trivial equilibrium solution where 


(s*L=° < 45 > 

which exists for all u/S and £l/S. Non-trivial equilibrium solutions where ( e/SK )„ ^ 0 
exist for intermediate ranges of tu/S and fi/S wherein the trivial equilibrium solution (45) 
typically becomes unstable. 

We will now show that the non-trivial equilibrium values of II^, IJ/oo, (633)00, and 
(7V e )°o are universal (i.e., completely independent of c o/S and f 1/S) for this hierarchy of 
models in plane homogeneous turbulent flows. A closed system of equations for the deter- 
mination of the temporal evolution of II, III, b 33 and e/SK are as follows: 


dll 

dt* 


dm 

dt* 


M-i)"- 


V e 2V e 

*° 33 — 7777 + 


eSK ZeSK 


^jK II + ^ In -^7jK 

V e 1 V £ 

+ ^7sk--^"7sk 


3e 

Jk 


(l - III + \u-4rz - II V 6 


2 e SK 


e SK 


-^7JK + 3 f' in jK + ^ ul jK 


+ ! A633 7 Tic - \ p ' u 77k 


+ I ftm ”7 Jk - W n 7 Jk 


(46) 


( 47 ) 
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dk 

dt* 


33 


£ 

~SK 


(>-?) 


'33 


IV e 
3 e SK 


+ 2^SK b33+ 2 ^ SK 3^0 


( 48 ) 


(49) 


Equations (46) - (47) are obtained by multiplying (37) with b and b 2 , respectively, and then 
taking the trace after the model (29) for % is substituted. The non-trivial equilibrium 
solutions are then obtained from the nonlinear algebraic equations 


+A«- + &m« - jA (|)^ + A(M» (j) 


(50) 


■I 1 - (f ) J (f ). - Ct). 

-(M» (f)_ + + ift/4 

+ i«‘“)-(7)--i W7 -(7)- 

+ iA/MM~(f)„-5ft^(|) oo = o 

t 1 “ (?) J (W “ “ 5 (f)oo + 5 A(M “ 

+5A [(ML - |n~] = 0 

(C.l - 1) (y) ^ - (<7,i - i) = 0 


(51) 


(52) 


(53) 


which are derived by setting the time derivatives on the left-hand-sides of (46) - (49) to 
zero and dividing by e/SK . The system of equations (50) - (53) will have solutions for 
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I loo, 1 1 loo, ( 633)00 and ( V jt)oo that are completely independent of a i/S and Q/ S - and hence 
universal - for plane homogeneous turbulent flows. 

This universal equilibrium structure of I loo, II loo, ( 633)00 and (P/e) 00 will now be utilized 
to obtain the simplest generic form of (29) that has the same equilibrium structure as the 
general model in the phase space of plane homogeneous turbulent flows. Due to these four 
universal invariants, the quadratic terms in the rapid pressure strain correlation are not 
linearly independent for plane homogeneous turbulent flows. This quadratic part Iljj' of (29) 
is defined as follows: 


= (3 5 K(bikbia5jl + bjkbtuSil — + @7 K{bik bktW , + bjkbtaWu)- (54) 

For plane homogeneous turbulent flows, a straightforward, although somewhat tedious, cal- 
culation yields the relationships 

2 III 

3 633 

ktbuW,, + b it b u W u = + Ij.TTa) (56) 

where we have made use of (38) and the fact that the anistropy tensor is of the form 


3* (55) 


2 2 
bikbkiSjl + bjkbklSxl — -bklbtm^mkSij = —633 (bik^jk + bjk^ik — 




( 611 &12 0 

&12 &22 0 

^ 0 0 &33 


(57) 


in such flows. Due to (55) - (56), and the fact that IIqo, IUoo, ( 633)00 and (P/£)oo are 
universal invariants for all plane homogeneous turbulent flows, it follows that the quadratic 
terms in the rapid pressure-strain are directly related to the linear terms in such flows. 
Consequently, the equilibrium structure of (29) in plane homogeneous turbulent flows will 
be indistinguishable from that of the substantially simpler model 

Tlij = c x ebij + c 2 e(bikb k j - - 6 mn 6 mn 6 ,y) 


+c 3 K!)ij + c 4 K(bik!>jk + bj k H>ik 


(58) 


— + c sK{b,kW jk + bjkWik) 

where Ci,c 2 , . . . ,c 5 are dimensionless constants and we have made use of the fact that 

tr(b-5) = -i| (59) 

tr-fb 1 S) = ii. 33 ^ (60) 
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which was also used in the derivation of (46) - (49). In alternative terms, (58) is topolog- 
ically equivalent to the general model (29) in so far as the equilibrium structure of plane 
homogeneous turbulent flows is concerned. 

It is rather striking that an analysis of the equilibrium states of arbitrary plane ho- 
mogeneous turbulence - coupled with the crucial physical constraint of invariance under 
coordinate transformations — collapses the general pressure-strain model 

n<,- = eAi(b) + KM ijU ( b)^, (61) 

(which can have as many as forty-five independent functions of b) to the substantially sim- 
plified model (58) that has only five undetermined constants. In the next section, a new 
model for the pressure-strain correlation will be developed. 

4. THE SSG MODEL: ITS ASYMPTOTIC ANALYSIS AND CALIBRATION 

Now, a new model for the pressure-strain correlation - which we will call the SSG model 
- will be developed based on the previous invariant dynamical systems analysis coupled with 
the following additional constraints: 

(i) Asymptotic consistency in the limit of small anisotropies 

(ii) Consistency with Rapid Distortion Theory for homogeneously strained turbulent flows 
that are initially isotropic 

(iii) Consistency with the equilibrium values for homogeneous shear flow obtained from the 
physical experiments of Tavoularis and Corrsin (1981) 

(iv) Consistency with the RDT results of Bertoglio (1982), for rotating shear flows, which 
predict that the most unstable flow occurs when the ratio of the rotation rate to the 
shear rate Q/S = 0.25 and that a flow restabilization occurs when Q/S > 0.5 

(v) Consistency with the results of physical experiments on the decay of isotropic turbu- 
lence and the return to isotropy of an initially anisotropic, homogeneous turbulence. 

Since the magnitude of the anisotropy tensor is relatively small (||b|| < 0.20 for most turbu- 
lent flows of engineering and scientific interest), we feel that terms which are of a comparable 
order of magnitude in bij should be maintained unless there is some overriding physical rea- 
son not to do so. In this fashion, the model can then be thought of as an asymptotically 
consistent truncation of a Taylor series expansion of Aij and Mijki in the variable 6~. Since 
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the simplified model for the rapid pressure-strain correlation in (58) is of 0(b), this suggests 
that c 3 - which in general can be a function of the invariants of b - should be replaced by 

o 3 - c;iii 

(where C 3 and are constants) for asymptotic consistency. Furthermore, since the model 
for the slow pressure-strain correlation is of 0 (b 2 ), and since most engineering flows have 
significant regions where "P > e, we will replace the constant Cj with the coefficient 

-(<?i + C{V/e) 


where C\ and C* are constants. This yields the following model for the pres sure- strain 


correlation: 


n»j = —(C\6 + C'V )6,j + Ciefiikbkj 


+ (C 3 - CjnijJfJi, 

_ 2 

-|- C4 K ( 6 t /c S jk + bjk~S%k X bmn^mn ) 


(62) 


+C,K(b, k W, l + bjtW*). 

Although (62) is topologically equivalent to (58) in so far as the equilibrium states are 
concerned, it will give rise to different temporal evolutions. We feel that it is better to 
use (62) as our final model for the pressure-strain correlation since all terms that are of a 
comparable order in bij have been maintained for asymptotic consistency. 

Before using constraints (ii) - (v) to calibrate the SSG model given by (62), a few com- 
ments are in order concerning the relationship between this new model and previously pro- 
posed models. The SSG model is not significantly more complicated than the commonly 
used Launder, Reece, and Rodi model which can be obtained from (62) in the linear limit as 
c*, C 2 and C 3 go to zero. In fact, the SSG model is substantially simpler than the recently 
proposed nonlinear models of Shih and Lumley (1985) and Fu, Launder, and Tselepidakis 
(1987) (see Appendix B). 

The coefficients C\ and C 2 have been calibrated by considering the return to isotropy 
problem (see Sarkar and Speziale 1989). Of course, for the return to isotropy problem, only 
the terms containing the coefficients Ci and C 2 in the pressure-strain correlation survive (i.e., 
the rapid pressure-strain correlation vanishes). Based upon realizability, dynamical systems 
considerations, and the phase space portrait of return to isotropy experiments, the following 
values of C\ and C 2 were arrived at by Sarkar and Speziale (1989): 


C x = 3.4 


(63) 
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C 2 = 3(Ci - 2) = 4.2. (64) 

Interestingly enough, the value of C\ — 3.4 is quite close to the value of 3.6 for the Rotta 
coefficient that is currently used in the basic model of Launder and his co-workers. However, 
as demonstrated by Sarkar and Speziale (1989), the quadratic term containing C 2 is crucial 
to properly capture the experimental trends. In Figure 1, the predictions of the SSG model 
and the Launder, Reece, and Rodi (LRR) model are compared in the £ — 77 phase space 
with the experimental data of Choi and Lumley (1984) for the return to isotropy from plane 
strain. The SSG model exhibits a curved trajectory that is well within the range of the 
experimental data; the LRR model - as well as any model for which C 2 = 0 - erroneously 
predicts a straight line trajectory. In Figures 2(a) - (b), the predictions of the SSG model 
and the LRR model for the temporal evolution of the anisotropy tensor are compared with 
experimental data for the relaxation from plane strain experiment of Choi and Lumley (1984) 
and plane contraction experiment of Le Penven, Gence, and Comte- Bellot (1985). The SSG 
model, on balance, yields improved predictions over the LRR model. For more detailed 
discussions and comparisons, the reader is referred to the paper by Sarkar and Speziale 
(1989) where this quadratic model for the slow pressure-strain correlation was compared 
with data from four independent experiments on the return to isotropy. 

Constraint (ii), which requires consistency with RDT for a homogeneously strained tur- 
bulence that is initially isotropic, is commonly enforced in the turbulence modeling literature. 
While the dynamical systems analysis presented in Section 3 can guarantee proper long-time 
behavior, this RDT constraint can be of considerable assistance in ensuring proper short- 
time behavior; if a model properly captures both the short and long-time behavior, it stands 
an excellent chance of performing well for all times. This RDT result requires that (see Crow 
1968) 

Hm % = (65) 

and, hence, that 

C3 = — (66) 

for the SSG model. We found that models which deviated significantly from (66) performed 
poorly in homogeneous shear flows (e.g., such models yielded spurious points of inflection in 
the time evolution of the turbulent kinetic energy). 

Constraints (iii) - (iv) were used to calibrate the remaining constants in the model, 
namely, C^, C^, C4, and C 5 as well as the constant C e 1 in the modeled e- transport equation. 
This was done using a value of C e 2 = 1.83 (as opposed to the more commonly adopted 
value of 1.92) which yields a power law decay in isotropic turbulence with an exponent of 
1.2 - a value which is in better agreement with available experimental data as discussed by 
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Reynolds (1987). It was not possible to obtain the exact equilibrium values of Tavoularis 
and Corrsin (1981) for homogeneous shear flow and satisfy the RDT results of Bertoglio for 
rotating shear flow (i.e., constraints (iii) - (iv)). A numerical optimization yielded the values 
of C* = 1.80, C7* = 1 . 30 ,C 4 = 1.25, C s = 0.40 and C eX = 1.44 as the best compromise. The 
equilibrium values obtained from the SSG model (using these values of the constants) for 
homogeneous shear flow are compared with the values obtained from the Launder, Reece, 
and Rodi model and the experiments of Tavoularis and Corrsin (1981) in Table 1. From 
these results, it is clear that the predictions of the SSG model are well within the range 
of the experiments whereas the predictions of the Launder, Reece, and Rodi model deviate 
significantly. Furthermore, the SSG model predicts that the largest growth rate in rotating 
shear flow occurs when fi/5 ~ 0.22 and that a flow restabilization occurs when fi / S > 0.53 
in comparison to the corresponding RDT results of fi / S — 0.25 and fi / S > 0.5 (see Bertoglio 
1982). These predictions of the SSG model are considerably better than those of the LRR 
model which erroneously predicts that the largest growth rate occurs when 9,/S = 0.14 
and that a flow restabilization occurs when O j S > 0.37. A more detailed discussion of the 
performance of the models in rotating shear flow will be presented in the next section. 

The SSG model has been carefully calibrated to perform well in shear flows both with 
and without added rotational strains. It is our belief that this will significantly enhance 
the performance of the model in turbulent boundary layers with streamline curvature - an 
analogous flow with a variety of important engineering applications. However, unlike other 
recently derived models for the pressure-strain correlation such as the Shih-Lumley (1985) 
model and the Fu, Launder, and Tselepidakis (1987) model, the SSG model given by equation 
(62) does not satisfy the strong form of realizability. The strong form of realizability (see 
Lumley 1978) constitutes a sufficient condition to guarantee positive component energies 
in homogeneous turbulent flows. The SSG model only satisfies a weak form of realizability 
wherein the turbulent kinetic energy is guaranteed to be positive; this is a direct consequence 
of the form of the modeled e-transport equation (see Speziale 1989b). We decided to opt for 
the weaker form of realizability for two major reasons. First, if the turbulent kinetic energy 
is positive, realizability can only be violated by fairly large anisotropies, such that 

IN > 3 

(where || • || is the L 2 norm or maximum eigenvalue) which are outside of the expected domain 
of applicability of such idealized models. It must be kept in mind that, so long as the model 
yields a positive turbulent kinetic energy, it can be applied to a flow (it is primarily negative 
kinetic energies that are computationally fatal). Second, it has been our experience that 
models which satisfy the strong form of realizability become computationally “stiff” in flows 
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with large anisotropies. This results from the fact that the finite difference form of the 
modeled equations usually do not exactly satisfy realizability (see Speziale and Mac Giolla 
Mhuiris 1989a). No such problems are encountered by the weak form of realizability since it is 
satisfied exactly by most standard numerical formulations of the model. Finally, it should be 
mentioned that the SSG model was not forced to satisfy material frame indifference (MFI) in 
the limit of two-dimensional turbulence (Speziale 1983) which constitutes another extreme 
constraint that is a rigorous consequence of the Navier-Stokes equations. It has recently 
become apparent to us that when such extreme constraints as MFI and strong realizability 
(correct as they may be for general flows) are applied to highly idealized models , there is a 
strong possibility that the model will become overly biased so that it performs poorly in the 
more commonly encountered turbulent flows. 


5. PERFORMANCE OF THE SSG MODEL IN HOMOGENEOUS FLOWS 

The SSG model given by equation (62) will now be tested in four independent homo- 
geneous turbulent flows. For the purpose of clarity, we will summarize the values of the 
constants that were arrived at in the previous section: 


Cl = 3.4, C\ = 1.80, C 2 = 4.2 

(67) 

C 3 = 7, Cl = 1.30, C 4 = 1.25 
5 

(68) 

's = 0.40, C e i = 1-44, C e2 = 1.83. 

(69) 


The problem of homogeneous turbulent shear flow in a rotating frame will be considered 
first. For this case, the mean velocity gradients and the rotation rate of the reference frame 
are given by 

(70) 

(71) 



respectively, in matrix form. The initial conditions correspond to a state of isotropic turbu- 
lence where 

bn = 0, K = K 0} e = e 0 (72) 

at time t = 0. It was shown by Speziale and Mac Giolla Mhuiris (1989a) that the solution 
only depends on the initial conditions through the dimensionless parameter €q/SKo ; the 
dependence of the solution on the shear rate and rotation rate is exclusively through the 
dimensionless parameter Q/*S\ Two types of equilibrium solutions have been established 
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for this problem (Speziale and Mac Giolla Mhuiris 1989a): one where (e/SK)^ = 0 which 
exists for all fi / S and one where ( ejSK)^ > 0 which exists only for a small intermediate 
band of values for fi/ S. The zero equilibrium value is associated predominantly with stable 
flow wherein K and e undergo a power-law time decay; the nonzero equilibrium values are 
associated with unstable flow wherein K and e undergo an exponential time growth. The two 
solutions undergo an exchange of stabilities for intermediate values of £l/S (which includes 
the case of pure shear where Q/S = 0). In this fashion, the second-order closures are able 
to account for both the shear instability - with its exponential time growth of disturbance 
kinetic energy - and the stabilizing (or destabilizing) effect of rotations on shear flow. 

In Figures 3(a) - (c), the predictions of the SSG model for the time evolution of turbulent 
kinetic energy are compared with those of the LRR model and the results of the large-eddy 
simulations of Bardina, Ferziger, and Reynolds (1983) for rotating shear flow. From these 
figures, it is clear that the SSG model does a much better overall job of capturing the trends 
of the large-eddy simulations. Several observations are noteworthy: (a) the LRR model 
exhibits too strong a growth rate for pure shear (£l/S = 0) in comparison to the SSG model 
and large-eddy simulations, (b) both the SSG and LRR models yield too weak a growth 
rate for the Q/S = 0.25 case, however, the SSG model is substantially better, and (c) the 
SSG model properly captures the weak growth rate that occurs for Q/S = 0.5, whereas the 
LRR model erroneously predicts a flow restabilization. The premature flow restabilization 
predicted by the LRR model at fl/S ~ 0.37 is somewhat serious since, in addition to the 
results of large-eddy simulations, linear stability theory and RDT predict that there should 
be unstable flow for the entire range of 0 < fI/5 < 0.5 (see Lezius and Johnston 1976 and 
Bertoglio 1982). As mentioned earlier, the SSG model does not predict a flow restabilization 
until a/S > 0.53. 

It would be useful at this point to compare the performance of the SSG model in rotating 
shear flow with that of some newer models that have been recently proposed. Three such 
models - the model of Shih and Lumley (1985), the model of Fu, Launder, and Tselepidakis 
(1987), and the RNG model of Yakhot and Orszag (1988) - were compared in a recent study 
of Speziale, Gatski, and Mac Giolla Mhuiris (1989). It was established in that study that 
the Fu, Launder, and Tselepidakis (FLT) model performed the best among these models 
in rotating shear flow. Hence, for simplicity, we will only compare the SSG model with 
the FLT model. In Figures 4(a) - (c), the results for the turbulent kinetic energy obtained 
from the SSG model and the FLT model for the rotation rates of jS = 0 ,ti/S = 0.25, 
and Q/S = 0.50 are compared with the large-eddy simulations of Bardina, Ferziger, and 
Reynolds (1983) for rotating shear flow. It is clear from these results that the SSG model 
properly captures the trends of the large-eddy simulations which indicate that all three cases 
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are unstable and that the Cl/S = 0.25 case has the strongest growth rate. On the otherhand, 
the FLT model erroneously predicts that the Cl/S = 0 and Cl/S = 0.25 cases are equally 
energetic and that the Cl/S = 0.5 case has undergone a restabilization. Comparable to the 
LRR model, the FLT model erroneously predicts a premature restabilization at Cl/S ~ 0.39. 
It may be of concern that a heavy emphasis has been placed on comparisons with large- 
eddy simulations for rotating shear flow (unfortunately, no direct simulations or physical 
experiments have been conducted for this problem). However, it must be emphasized that 
the critical evaluations have been based on which states should be more energetic - results 
which have been confirmed independently by RDT and linear stability theory. 

A bifurcation diagram for the general hierarchy of closure models (61) is shown in Figure 
5 for rotating shear flow. Here, the equilibrium value of (e/SK) oa is plotted as a function of 
Cl/S. The SSG model as well as the other commonly used models have the same topological 
structure in rotating shear flow as indicated in Figure 5. There axe two equilibrium solutions: 
the solution where (e/ SK) «, = 0 exists for all Cl/S but becomes unstable in the interval AB\ 
the nonzero solution for (e/5iif) 00 , which lies on the semi-ellipse ACB, exchanges stabilities 
with the trivial solution ( e/SK )<» = 0 in the interval A < Cl/S < B. For Cl/S < A — SA and 
Cl/S > B -f SB (where 6 A and SB represent a small increment that depends on the model) 
the trivial equilibrium value of (e/SK)^ = 0 is associated with solutions where the kinetic 
energy undergoes a power law decay with time; for A — SA < Cl/S < B + SB, this trivial 
solution is associated with solutions where the kinetic energy undergoes a power law growth 
with time. The nonzero equilibrium values (e/SK) „ > 0 (on the semi-ellipse ACB ) are 
associated with solutions where the kinetic energy grows exponentially with time. It can be 
shown (see Speziale and Mac Giolla Mhuiris 1989a) that the growth rate A for A < Cl/S < B 
is given by 

x (73) 

where a = (C e 2 — l)/(C e i — 1). Hence, point C - which corresponds to the maximum value 
of ( e/SK)oo - represents the most energetic state with the largest growth rate of kinetic 
energy. 

The coordinates [Cl/ S,(e/ SK) <*,] of points A, B, and C (in Figure 5) for the Launder, 
Reece, and Rodi model and the SSG model are given below: 


LRR Model 


A = [-0.09,0], B = [0.37,0], C = [0.14,0.167] 


SSG Model 


A = [-0.09,0], B = [0.53,0], C = [0.22,0.254]. 
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The improved performance of the SSG model in rotating shear flow is largely due to the 
fact that its most energetic state (point C on the bifurcation diagram shown in Figure 5) is 
located in close proximity to Cl/S = 0.25 - the value predicted by rapid distortion theory. 
However, it needs to be mentioned at this point that the reason we were not able to satisfy 
this RDT result exactly is due to a defect in the general hierarchy of models (61). Due to 
(73) and the fact that the bifurcation diagram is symmetric about its most energetic state 
(point C in Figure 5), the general hierarchy of models (61) erroneously predicts Richardson 
number similarity if point C is located at Cl/S = 0.25. Such models will yield solutions for 
K and e that scale with the Richardson number 


Ri = -2(fi/5)(l - 2fi/5) 


(74) 


and, thus, erroneously predict that the Cl/S = 0 and Cl/S = 0.5 cases are identical. Large- 
eddy simulations, RDT (Bertoglio 1982), and independent mathematical analyses of the 
Navier-Stokes equations (Speziale and Mac Giolla Mhuiris 1989b) indicate that the Cl/S = 0 
and Cl/S = 0.5 cases are distinct. By moving the most energetic state a small distance to 
the left of Cl/S — 0.25 - as is done with the SSG model - the proper growth rates obtained 
from large-eddy simulations for Cl/S = 0 and Cl/S = 0.5 can be reproduced (see Figures 3(a) 
and 3(c)). However, the substantially larger growth rate for Cl/S = 0.25 shown in Figure 
3(b) (which has independent support in the RDT calculations of Bertoglio 1982), cannot be 
reproduced by the SSG model. This is a defect in the general hierarchy of models (61) that is 
intimately tied to their prediction of universal equilibrium values for //<» , ///qo , ( 633)00 and 
('P/e) 00 in plane homogeneous turbulent flows - an oversimplification that is not supported 
by physical or numerical experiments. 1 Nonetheless, despite this deficiency, the SSG model 
performs reasonably well - and is superior to other existing second-order closures - for 
rotating shear flow as evidenced by Table 1, Figures 3(a) - (c), and Figures 4(a) - (c). 

Now, we will examine the performance of the SSG model in homogeneous plane strain 
turbulence for which the mean velocity gradients are given by 


dvi 
dx j 


S 0 0 \ 

0-5 0 

0 0 0 / 


(75) 


and the turbulence evolves from an initial state of isotropy. Comparisons of the model predic- 
tions will be made with the direct numerical simulations of Lee and Reynolds (1985) on plane 
strain. Such comparisons must be made with caution due to the low turbulence Reynolds 
numbers of the direct simulations. However, comparisons with physical experiments (e.g., 

*It is not possible to tie this deficiency to the modeled e - transport equation since all dependence on e 
can be eliminated in the RDT limit. 
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Tucker and Reynolds 1968 and Gence and Mathieu 1980) are equally problematical due to 
the uncertainty in the initial conditions for e/SK and possible large-scale contamination 
from the walls of the test apparatus. 

In Figure 6, the time evolution of the turbulent kinetic energy for the LRR model and SSG 
model are compared with the direct simulations of Lee and Reynolds (1985) for plane strain 
corresponding to the initial condition eo/SKo = 2.0. It is clear from this figure that both 
models perform extremely well. In Figure 7, the time evolution of the non-zero components 
of the anisotropy tensor are shown. Although the quantitative accuracy of the models is 
not extremely good, it is clear that the SSG model does better than the LRR model and 
reproduces the crucial trends of the direct simulations. In Figures 8-9, the time evolution 
of the turbulent kinetic energy and non-zero components of the anisotropy tensor axe shown 
corresponding to the initial condition of eo/SKo = 1.0. The same conclusions can be drawn 
from these results: the SSG model yields improved predictions over the LRR model and, 
on balance, compares reasonably well with the direct simulations which would be expected 
to have somewhat elevated anisotropies due to the lower turbulence Reynolds number. We 
will not make more extensive comparisons with the predictions of other turbulence models 
since our main purpose here was to simply establish that the alterations made in the LRR 
model - to yield the SSG model with its improved behavior in rotating shear flows - do not 
compromise its performance in plane strain. 

Finally, we will examine the performance of the SSG model for the axisymmetric contrac- 
tion and expansion in homogeneous turbulence. Since the SSG model (like virtually all other 
existing models for the pressure-strain correlation) was calibrated based on plane homoge- 
neous turbulent flows, it would be desirable to assess its performance in a three-dimensional 
flow. For the axisymmetric contraction, the mean velocity gradients are given by 


dvi 

dxj 


S 

0 

0 


-is 

0 


0 

0 

-is 


whereas in the axisymmetric expansion they take the form 


dvi 

dxj 


“J- 


(76) 


(77) 


The time evolution of each of these turbulent flows - from an initially isotropic state - will 
be considered. Hence, as with plane shear and plane strain, the solutions will only depend 
on the initial conditions through the parameter e 0 / SK 0 . Comparisons will be made with 
the predictions of the LRR model and the direct numerical simulations of Lee and Reynolds 
(1985) for the same reasons as cited above. 
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In Figure 10, the time evolution of the turbulent kinetic energy for the axisymmetric 
contraction is shown corresponding to the initial condition Eo/SKq = 0.179 which was taken 
from the direct simulations of Lee and Reynolds (1985). From these results, it is clear that 
the SSG model yields noticeably improved predictions over the LRR model; however, both 
models predict growth rates that are smaller than those in the direct simulations. In Figure 
11, the time evolution of the nonzero components of the anisotropy tensor are Bhown for the 
axisymmetric contraction where Eq/SKq = 0.179. While the differences between the SSG 
and LRR models is small, it is clear that the SSG model yields results that are more in line 
with the direct simulations. 

In Figure 12, the time evolution of the turbulent kinetic energy for the axisymmetric 
expansion is shown for the initial condition Eq/SKq = 2.45 taken from the direct numerical 
simulations of Lee and Reynolds (1985). It is clear from this figure that both the SSG 
and LRR models yield results that are in excellent agreement with the direct simulations. 
However, the time evolution of the nonzero components of the anisotropy tensor shown in 
Figure 13 show more significant differences. Here, the predictions of the SSG model appear 
to be substantially better than those of the LRR model. 

6. CONCLUSIONS 

The modeling of the pressure-strain correlation of turbulence has been considered based 
on invariance arguments and a dynamical systems approach. Several important conclusions 
have been drawn about the standard hierarchy of closures (61) which led to the development 
of a new model - the SSG model. A summary of these findings can be given as follows: 

(i) It was proven that the standard hierarchy of models yields non-trivial values for the 
equilibrium states 11^, III^, (633)00, and ( V/e)^ that are universal (i.e., that do not 
depend on u/S, ti/S or the initial conditions) for plane homogeneous turbulent flows. 
As a direct consequence of these universal invariants, it was shown that, for plane 
homogeneous turbulent flows, the general model (61) for the pressure-strain correlation 
is topologically equivalent to a substantially simpler model — the SSG model - which 
is only quadratically nonlinear in the anisotropy tensor. 

(ii) The SSG model was calibrated by using existing data from isotropic decay experiments, 
return to isotropy experiments, and homogeneous shear flow experiments along with 
the RDT results of Crow (1968) and Bertoglio (1982). By means of this more systematic 
method of calibration, the SSG model was demonstrated to perform better than the 
Launder, Reece, and Rodi model - as well as the newer models of Shih and Lumley 
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and Fu, Launder, and Tselepidakis - for a variety of homogeneous turbulent flows. 
The flows that were examined included the challenging test case of rotating shear flow 
(where rotations can either stabilize or destabilize the flow) and the axisymmetric 
expansion/contraction which constitutes a three-dimensional mean turbulent flow. 

(iii) Although the SSG model performs reasonably well for a variety of homogeneous tur- 
bulent flows, there are still major deficiencies with it that are intrinsic to this general 
hierarchy of models. These deficiencies, emanate from the prediction of universal equi- 
librium values for /ioo, ///«,, ( 633)00 and ('P/e) 0O in plane homogeneous turbulent flows 
- an obvious oversimplification that is not supported by physical experiments. As a 
result of this deficiency, the general model (61) erroneously predicts that rotating shear 
flow has growth rates that are symmetric about their most energetic value. Hence, in 
order to satisfy the RDT constraint of Bertoglio (1982) — which puts the most energetic 
state at f 1/S = 0.25 - the models must exhibit Richardson number similarity. This is 
inconsistent with the Navier-Stokes equations as proven by Speziale and Mac Giolla 
Mhuiris (1989b) and illustrated by Bardina, Ferziger, and Reynolds (1983). 

(iv) Since the general model (61) for the pressure-strain correlation gives an incomplete 

picture of plane homogeneous turbulent flows no matter what form is taken for Aj(b) 
and we feel that the process of adding highly complex nonlinear terms in b XJ 

is somewhat questionable. Such complex nonlinear terms in the rapid pressure-strain 
correlation have been largely motivated by the desire to satisfy the strong form of 
realizability. However, it must be remembered that the strong form of this constraint 
only constitutes a sufficient condition for the satisfaction of realizability in homoge- 
neous turbulent flows. Due to the relatively simple topological structure of the general 
model in rotating shear flow - which is in no way altered by the addition of more 
complex nonlinearities in K - the application of the strong form of realizability either 
removes the degrees of freedom necessary to properly calibrate the model or leads to 
stiff behavior. 

Despite the deficiencies discussed above, the SSG model seems to perform moderately 
well in a variety of homogeneous turbulent flows as documented in this paper. While further 
improvements would be most welcome, we feel that it is unlikely that they will come from the 
standard hierarchy of models (61). Fundamentally new approaches will be needed. Future 
research will be directed on two fronts. The SSG model will be implemented in a full 
second-order closure for the computation of a variety of complex aerodynamic flows that are 
of technological interest. We believe that when the SSG model is used within the framework 
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of a sound second-order closure, it may be possible to obtain acceptable engineering answers 
for a range of turbulent shear flows with streamline curvature. In parallel with this effort, 
entirely new directions in modeling the pressure-strain correlation will be pursued. These 
will involve the introduction of a tensor length scale - to better account for anisotropies - 
and the possible solution of a transport equation for M+jki to account for history effects in 
the rapid pressure- strain. A closer examination of these issues will be the subject of a future 
paper. 
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APPENDIX A 


Consider the tensor function 

= Fiitb.S',*'). (Al) 

Form invariance under coordinate transformations (28) requires that be of the form (see 
Smith 1971): 

F = Qol + c^b + 02b 2 + Q3S + 0:48 

+a 5 (bS' + S'b) + a 6 (b 2 S' + S'b 2 ) 

+a 7 (bS' 2 + S ,2 b) + a 8 (b 2 S' 2 + S% 2 ) 

+a 9 (bw / — uT'b) + a^o^bur* ( A2) 

+au(b 2 aJ v — ZT'b 2 ) + ai2(w" , buT' 2 — uT ,2 bw J ) 

+a 13 (S — uT^S ) + ax4oT / S U 1 

+ai 5 (S 2 uT / — Z0S 2 ) + ai6(t*r / S U 42 — uT /2 S Z?) + ocnZT 12 
where depend on all possible invariants, i.e., 

ax = «<(//, III) tr ZT 12 , tr S >2 ,tr S' 3 ,fr bS ',tr b 2 S', 

tr bS' 2 ,tr b 2 S' 2 ,tr bZ0 2 >tr b 2 ^ 2 ,^ ZfbZJ 12 , 

(A3) 

tr SV 2 , tr ZT'S'ZZ 12 , tr bS'aT'jir b 2 S'aT', 

tr b^ 2 S', tr bS'V), % = 0, 1, 2 , .... 17. 

Taking the linear part of Fa in S[j and ZJ'- yields the equation 

F< L > = /3 0 I + Ab + /3 2 b 2 +/3 3 S' + /3 4 (bS' 

+S'b) + /3 5 (b 2 S' + S'b 2 ) (A4) 

+0 6 (bZ0 — uT'b) + 0 7 (b 2 Z0 — u^b 2 ) 

where 

00 = Poo (II, III) + 0oi(H, III)tr( b • S') + 002 (II, III)tr( b 2 • s') (A5) 

01 = 0io(H, HI) + 0u(II, III)tr(b • S') + 0 12 (II, III)tr(b 2 • s') (A6) 
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fa = fao(II, III ) + Mil, III)tr(b ■ S') + fa 2 {II, III)tr{ b 2 • S') 
fa = fa(II,III), * = 3,4, • • • ,7. 


(A7) 

(A8) 


Then, by taking 
obtain II jj : 


the deviatoric part of (since Ilij is traceless) and multiplying by e we 



(A9) 


Equation (29) results when (A4) - (A8) are substituted into (A9) and the identities in 
equation (27) are made use of. 


27 


APPENDIX B 


The pressure-strain models of Shih and Lumley (1985) and Fu, Launder, and Tselepidakis 
(1987) take the following form: 

Shih-Lumley Model 


where 


n* = 


3 (! a 5 ) jr* ^ +e mkj n 


—fitbij + ^- + 8a 5 ^ K~5ij 

B 

+ (I + f° 5 ) h + 


(m % 

+Tjk + emkiQ 


Tmn&mn &ij 


+Tjk [dx] + 3 


TtyitV ^ mn j 


6 4 

“ 1 “ y^T \^~ik jk * 4 * 

5 k { T,fc {eft + em ^ n "*) 


+ T jk ^ T kl 2T^Tj' m 3 ? ^ m j> 


a s = ~(1 + O.sVf), F = 1 + 911+ 27/77 
77 = -^,%, 777 = ^bijbjkbid 

(3 = 2 + ip{80.Un[l + 62.4(— 77 + 2.3777)]} 


( 51 ) 


( 52 ) 

( 53 ) 

( 54 ) 
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+r{l6II{T ik W jk + TjkWik) 

24 b m ib n j(T m )cW n k 4 " 7~nk mk ) }" 

where 

C* = —6011, 7 = 0.6, r = 0.7. (£6) 
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Equilibrium 

Values 

LRR Model 



(^ll)oo 

0.193 

0.204 

0.201 

(^ 22)00 

-0.096 

-0.148 

-0.147 

(^ 12)00 

-0.185 

-0.156 

-0.150 

(SK/e) n 

5.65 

5.98 

6.08 


Table 1. Comparison of the predictions of the Launder, Reece and 
Rodi (LRR) model and the Speziale, Sarkar and Gatski (SSG) model 
with the experiments of Tavoularis and Corrsin (1981) on homogene- 
ous turbulent shear flow. 
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Phase space of the return to isotropy problem: Comparison of the predic- 
tions of the SSG model and LRR model with the plane strain experiment 

of Choi and Lumley (1984). 
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Figure 2. 


Time evolution of the anisotropy tensor during the return to isotropy. Com- 
parison of the predictions of the SSG model and LRR model with experi- 
ments: (a) the plane strain experiment of Choi and Lumley (1984), and (b) 
the plane contraction experiment of Le Penven, Gence, and Comte-Bellot 
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Figure 2. Time evolution of the anisotropy tensor during the return to isotropy. Com- 

parison of the predictions of the SSG model and LRR model with experi- 
ments: (a) the plane strain experiment of Choi and Lumley (1984), and (b) 
the plane contraction experiment of Le Penven, Gence, and Comte-Bellot 
(1985). 
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Figure 3. Time evolution of the turbulent kinetic energy in rotating shear flows. 

Comparison of the predictions of the SSG model and LRR model with 
the large-eddy simulation of Bardina, Ferziger, and Reynolds (1983) for 
eo/SKo = 0.296: (a) Cl/S = 0, (b) Cl/S = 0.25, and (c) Cl/S = 0.5. 
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Figure 3. Time evolution of the turbulent kinetic energy in rotating shear flows. 

Comparison of the predictions of the SSG model and LRR model with 
the large-eddy simulation of Bardina, Ferziger, and Reynolds (1983) for 
e 0 /SK 0 = 0.296: (a) f 'l/S = 0, (b) fi/S = 0.25, and (c) Q/S = 0.5. 
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LRR Model 


SSG Model 



Time evolution of the turbulent kinetic energy in rotating shear flows. 
Comparison of the predictions of the SSG model and LRR model with 
the large-eddy simulation of Bardina, Ferziger, and Reynolds (1983) for 
eo/SKo = 0.296: (a) f 1/S = 0, (b) fl/S = 0.25, and (c) Sl/S = 0.5. 
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Comparison of the models with large-eddy simulations (LES) of rotating 
shear flow for Q/S = 0,0.25 and 0.5: (a) LES of Bardina, Ferziger, and 
Reynolds (1983), (b) FLT model, and (c) SSG model. 
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Comparison of the models with large-eddy simulations (LES) of rotating 
shear flow for Cl/S = 0,0.25 and 0.5: (a) LES of Bardina, Ferziger, and 
Reynolds (1983), (b) FLT model, and (c) SSG model. 
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Comparison of the models with large-eddy simulations (LES) of rotating 
shear flow for £l/S = 0,0.25 and 0.5: (a) LES of Bardma, Ferziger, and 
Reynolds (1983), (b) FLT model, and (c) SSG model. 
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PLANE STRAIN 


Time evolution of the turbulent kinetic energy in plane strain for e 0 /SK 0 - 
2.0. Comparison of the predictions of the LRR model and SSG model with 
the direct simulations of Lee and Reynolds (1985). 
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Figure 7. Time evolution of the anisotropy tensor in plane strain for £ 0 /SK 0 = 2.0. 

Comparison of the predictions of the LRR model and SSG model with the 
direct simulations of Lee and Reynolds (1985). 
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Figure 8. Time evolution of the turbulent kinetic energy in plane strain for e 0 /SK 0 = 

1.0. Comparison of the predictions of the LRR model and SSG model with 
the direct simulations of Lee and Reynolds (1985). 
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Figure 9. Time evolution of the anisotropy tensor in plane strain for £o/SKq - 1.0. 

Comparison of the predictions of the LRR. model and SSG model with the 
direct simulations of Lee and Reynolds (1985). 
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Figure 10. 


Time evolution of the turbulent kinetic energy in the axisymmetric con- 
traction for eo/SKo = 0.179. Comparison of the predictions of the LRR 
model and SSG model with the direct simulations of Lee and Reynolds 
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Figure 11. Time evolution of the anisotropy tensor in the axisymmetric contraction 

for eo/SKo = 0.179. Comparison of the predictions of the LRR model and 
SSG model with the direct simulations of Lee and Reynolds (1985). 
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Figure 12. 


Time evolution of the turbulent kinetic energy in the axisymmetric expan- 
sion for Eo/SKo = 2.45. Comparison of the predictions of the LRR model 
and SSG model with the direct simulations of Lee and Reynolds (1985). 
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Figure 13. Time evolution of the anisotropy tensor in the axisymmetric expansion for 

eo/SKo = 2.45. Comparison of the predictions of the LRR model and SSG 
model with the direct simulations of Lee and Reynolds (1985). 
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